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Abstract 



We consider the problem of approximating all real roots of a square-free poly- 
nomial /. Given isolating intervals, our algorithm refines each of them to a width 
of 2~ L or less, that is, each of the roots is approximated to L bits after the binary 
point. Our method provides a certified answer for arbitrary real polynomials, only 
considering finite approximations of the polynomial coefficients and choosing a 
suitable working precision adaptively. In this way, we get a correct algorithm that 
is simple to implement and practically efficient. Our algorithm uses the quadratic 
interval refinement method; we adapt that method to be able to cope with inaccu- 
racies when evaluating /, without sacrificing its quadratic convergence behavior. 
We prove a bound on the bit complexity of our algorithm in terms of degree, size 
and separation of the roots, that is, parameters exclusively related to the geometric 
location of the roots. Our bound improves previous work on integer polynomials 
by a factor of deg/ and essentially matches best known theoretical bounds on root 
approximation which are obtained by very sophisticated algorithms. 

1 Introduction 

The problem of computing the real roots of a polynomial in one variable is one of the 
best studied problems in mathematics. If one asks for a certified method that finds all 
roots, it is common to write the solutions as a set of disjoint isolating intervals, each 
containing exactly one root; for that reason, the term real root isolation is common in 
the literature. Simple, though efficient methods for this problem have been presented, 
for instance, based on Descartes' rule of signs [6], or on Sturm's theorem [7 1. Recently, 
the focus of research shifted to polynomials with real coefficients which are approx- 
imated during the algorithm. It is worth to remark that this approach does not just 
generalize the integer case but has also lead to practical iflOl [161 and theoretical 1(171 
improvements of it. 

We consider the related real root refinement problem: assuming that isolating inter- 
vals of a polynomial are known, refine them to a width of 2~ L or less (where L > is 
an additional input parameter). Clearly, the combination of root isolation and root re- 
finement, also called strong root isolation, yields a certified approximation of all roots 
of the polynomial to an absolute precision of 2~ L or, in other words, to L bits after the 
binary point in binary representation. 

We present a solution to the root refinement problem for arbitrary square-free poly- 
nomials with real coefficients. Most of the related approaches are formulated in the 
REAL-RAM model where exact operations on real numbers are assumed to be avail- 
able at unit costs. In contrast, our approach considers the coefficients as bitstreams, that 
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is, it only works with finite prefixes of its binary representation, and we also quantify 
how many bits are needed in the worst case. The refinement uses the quadratic interval 
refinement method [ 1 1 (QIR for short) which is a quadratically converging hybrid of the 
bisection and secant method. We adapt the method to work with an increasing working 
precisions and use interval arithmetic to validate the correctness of the outcome. In 
this way, we obtain an algorithm that always returns a correct root approximation, is 
simple to implement on an actual computer (given that arbitrary approximations of the 
coefficients are accessible), and is adaptive in the sense that it might succeed with a 
much lower working precision than asserted by the worst-case bound. 

We provide a bound on the bit complexity of our algorithm. To state it properly, 
we first define several magnitudes depending on the polynomial which remain fixed 
throughout the paper. Let 

d 

/(x) g R[x] (1.1) 

1=0 

be a square-free polynomial of degree d > 2 with \a ( /\ > 1 andT:= |~log(max,- |fl,-|)] > 1. 
We denote the roots of / by zi,---,Zd, and, w.l.o.g., we can assume that the roots are 
numbered such that the first m roots z\ , . . . ,z m are all the real roots of /. For each z,, 
Oi = a(zi,f ) := min^,- \zj — Zj\ denotes the separation of zu : = L"=i logo) -1 an ^ 
:= log(max,- |z,|) the logarithmic root bound of /. An interval / = (a,b) is called 
isolating for a root zt if / contains zt and no other root of F. We set mid(7) = for 
the center and w(I) := b — a for the width of I. 

Main Result. Given initial isolating intervals for the roots of f, our algorithm 
refines one interval to width 7r L using 

d(d(dT f + T.f) 2 +dL) 

bit operations and refines all intervals using 

d(d(dr f + Y. f ) 2 +d 2 L) 

bit operations, where O means that we ignore logarithmic factors. To do so, our algo- 
rithm requires the coefficients off in a precision of at most 

d(dT f +Lf+L) 

bits after the binary point. 

For the analysis, we divide the sequence of QIR steps in the refinement process 
into a linear sequence where the method behaves like bisection in the worst case, and 
a quadratic sequence where the interval is converging quadratically towards the root, 
following the approach in [ 1 1 1. We do not require any conditions on the initial intervals 
except that they are disjoint and cover all real roots of F; an initial normalization phase 
modifies the intervals to guarantee the efficiency of our refinement strategy. 

We remark that, using the recently presented root solver from ifTTIl . obtaining initial 
isolating intervals can be done with 0{d{dTf + £/) 2 ) bit operations using coefficient 
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approximations of / to 0(dYf +£/) bits after the binary point. Combined with the 
latter result on root isolation, our complexity result also gives a bound on the strong 
root isolation problem. 

The case of integer coefficients is often of special interest, and the problem has 
been investigated by previous work [11] for this restricted case. In the latter work, the 
complexity of root refinement was bounded by 0(d A x 2 +d 3 L). We improve this bound 
to 

6(dh 2 +d 2 L). 

The difference in the complexities is due to a different approach to evaluate the sign 
of / at rational points which is the main operation in the refinement procedure: for an 
interval of size 2~ l , the evaluation of / at the endpoints of the interval has a complexity 
of 0(d 2 (z + £)) when using exact rational arithmetic because evaluated function values 
can consist of up to d{x + £) bits. However, we show that we can still compute the sign 
of the function value with certified numerical methods using the substantially smaller 
working precision of 0(dt + £). We remark that the latter result certainly only applies 
to points whose distance to a root is not much smaller than 2 , thus, we modified the 
QIR method in way such that the latter requirement is given. 

Related work. The problem of accurate root approximation is omnipresent in math- 
ematical applications; certified methods are of particular importance in the context of 
computations with algebraic objects, for instance, when computing the topology of 
algebraic curves |j5]|9l or when solving systems of multivariate equations J2). 

The idea of combining bisection with a faster converging method to find roots of 
continuous functions has been first introduced in Dekker's method and elaborated in 
Brent's method; see [4| for a summary. However, these approaches assume exact arith- 
metic for their convergence results. 

For polynomial equations, numerous algorithms are available, for instance, the 
Jenkins-Traub algorithm or Durant-Kerner iteration; although they usually approxi- 
mate the real roots very fast in practice, general worst-case bounds on their arithmetic 
complexity are not available. In fact, for some variants, even termination cannot be 
guaranteed in theory; we refer to the survey [ 15 1 for extensive references on these and 
further methods. 

The theoretical complexity of root approximation has been investigated by Pan lfl4ll . 
Assuming all roots to be in the unit disc, he achieves a bit complexity of d(r? +n 2 L) 
for approximating all roots to an accuracy of 2~ L , which matches our bound if L is 
the dominant input parameter. His approach even works for polynomials with multiple 
roots. However, as Pan admits in 1 15 1, the algorithm is difficult to implement and so is 
the complexity analysis when taking rounding errors in intermediate steps into account. 
Moreover, it appears unclear whether his bound can be improved if only a single root 
needs to be approximated. 

We improve on the first version of this paper |[T2l in two ways: first of all, in our 
bit complexity result, we remove the dependence on the coefficient size and, thus, re- 
late the hardness of root approximation to parameters that exclusively depend on the 
geometric location of the roots; we shortly expose in Section[6]who to benefit from this 
approach. Also, in this work, we redefine the threshold for the interval width that guar- 
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Algorithm 1 Eqir: Exact Quadratic Interval Refinement 



Input: / e R[x] squai-e-free, / = (a,b) isolating, iV = 2 2 ' eN 
Output: (J,N r ) with J C / isolating for E, and W'eN 
l: procedure eqir(/,7 = (a,b),N) 

ifN = 2, return (Bisection(/,/),4). 
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rri <r- a + round(iV /(a ^ fl j w )t» > m' w a + j^l^(b - a) 

s <- sign(/(m')) 

if s = 0, return ([to',to'],°°) 

if 5 = sign(/(a)) and sign(/(m' + £»)) = sign (/(&)), return ([m^m' + ffl],^ 2 ) 
if 5 = sign(/(&)) and sign(/(m' — to)) = sign(/(a)), return Qm' — to,m'],N 2 ) 
Otherwise, return (I,\/N). 
end procedure 



antees quadratic convergence (Defintion 12 1; in this way, we get rid of the magnitude 
R = log |res(/,/')| , which is a pure artifact of the analysis of lfl2l . 

Outline. We summarize the (exact) QIR method in Section[2] A variant using only 
approximate coefficients is described in Section [3] Its precision demand is analyzed in 
Section |4] Based on that analysis of a single refinement step, the complexity bound of 
root refinement is derived in Section[5] We end with concluding remarks in Section|6] 

2 Review on exact QIR 

Abbott's QIR method HJQjD * s a hybrid of the simple (but inefficient) bisection method 
with a quadratically converging variant of the secant method. We refer to this method 
as EQIR, where "E" stands for "exact" in order to distinguish from the variant presented 
in Section [3] Given an isolating interval / = (a,b) for a real root | of /, we consider 



the secant through (a,f(a)) and (b, f(b)) (see also Figure 3.1 1. This secant intersects 
the real axis in the interval /, say at x-coordinate m. For / small enough, the secant 
should approximate the graph of the function above / quite well and, so, m f=s t, should 
hold. An Eqir step tries to exploit this fact: 

The isolating interval / is (conceptually) subdivided into jV subintervals of same 
size, using 1 equidistant grid points. Each subinterval has width CO := Then 
m', the closest grid point to m, is computed and the sign of f(m') is evaluated. If that 
sign equals the sign of /(a), the sign of f(m' + co) is evaluated. Otherwise, f(m! — to) is 
evaluated. If the sign changes between the two computed values, the interval (rri ,m' + 
to) or the interval (to' — 0),to'), respectively, is set as new isolating interval for In 
this case, the EQIR step is called successful. Otherwise, the isolating interval remains 
unchanged, and the Eqir step is called failing. See Algorithm [T] for a description in 
pseudo-code. 

In fin , the root refinement problem is analyzed using the just described EQIR 
method for the case of integer coefficients and exact arithmetic with rational numbers. 
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Figure 3.1: Illustration of an Aqir step for N = 4. 



For that, a sequence of Eqir steps is performed with N = 4- initially. After a successful 
Eqir step, N is squared for the next step; after a failing step, N is set to y/N. If N 
drops to 2, a bisection step is performed, and N is set to 4 for the next step. In ifTTIl . a 
bound on the size of an interval is provided to guarantee success of every Eqir and, 
thus, quadratic convergence of the overall method. 



3 Approximate QIR 

The most important numerical operation in an Eqir step is the computation of f(xo) 
for values xq E I. Note that f(xo) is needed for determining the closest grid point m' 
to the secant (Step [4] of Algorithm [TJ, and its sign is required for checking for sign 
changes in subintervals (Steps |6]|8]l. 



What are the problems if / is a bitstream polynomial as in (1.1 1, so that /(jcq) 



can only be evaluated up to a certain precision? First of all, j^yzjnp. can only be 
computed approximately, too, which might lead to checking the wrong subinterval in 
the algorithm if m is close to the center of a subinterval. Even more seriously, if f(xo) 
is zero, then, in general, its sign can never be evaluated using any precision. Even if 
we exclude this case, the evaluation of f(xo) can become costly if xq is too close to a 
root of /. The challenge is to modify the QIR method such that it can cope with the 
uncertainties in the evaluation of /, requires as few precision as possible in a refinement 
step and still shows a quadratic convergence behavior eventually. 

Bisection is a subroutine called in the QIR method if N = 2; before we discuss the 
general case, we first describe our variant of the bisection in the bitstream context. Note 
that we face the same problem: Writing mid(7) as the center of / = (a,b), /(mid(7)) 
might be equal or almost equal to zero. We will overcome this problem by evaluating 
/ at several x-coordinates "in parallel". For that, we subdivide / into 4 equally wide 
parts using the subdivision points mj := a + j ■ for 1 < j < 3. We also assume that 
the sign of / at a is already known. We choose a starting precision p and compute 
f{m\), . . . ,f{mo,) using interval arithmetic in precision p (cf. Section|4]for details). If 
less than 2 out of 3 signs have been determined using precision p, we set p 2p and 



Algorithm 2 Approximate Bisection 



Input: / e R[x] square-free, / = (a,b) isolating, s = sign(/(a)) 
Output: J CI isolating with 2 • w(J) < w(I). 

i: procedure Approximate_Bisection(/,/= {a,b),s) 

2: V^[a+(i-l)-^,i = l,...,5] 

3: S= [s,0,0,0,-s] 

4: p 4-2 

5: while 5 contains more than one zero do 
6: for i=2,. . . ,4 do 

7: If S\i] = 0, set S[i] 4- sign«8(/(V[i]),p) 

8: end for 

9: P<-2p 

10: end while 

11: Find v,w, such that S[v] ■ S[w] =-lA (v+ 1 = vvV (v + 2 =wA5[v+l] =0)) 
12: return (V[v],V[w]) 
13: end procedure 



repeat the calculation with increased precision. Once the sign at at least 2 subdivision 
points is determined, we can determine a subinterval of at most half the size of / that 
contains ^ (Algorithm [2]). We will refer to this algorithm as "bisection", although the 
resulting interval can also be only a quarter of the original size. Note that / can only 
become zero at one of the subdivision points which guarantees termination also in the 
bitstream context. Moreover, at least 2 of the 3 subdivision points have a distance of 
at least to This asserts that the function value at these subdivision points is 
reasonable large and leads to an upper bound of the required precision (Lemma[5]). 

We next describe our bitstream variant of the QIR method that we call approximate 
quadratic interval refinement, or AQIR for short (see also Figure |3~T| for the illustration 
of an Aqir step for N = 4). Compared to the exact variant, we replace two substeps. 
In Step |4j we replace the computation of X :— ^ /(Jyff^b) as f°ll° ws: F° r a working 
precision p, we evaluate f(a) and f{b) via interval arithmetic with precision p (blue 
vertical intervals in the above figure) and evaluate jV {^"L, > with interval arithmetic 



accordingly (cf. Section|4|). Let J = (c,d) denote the resulting interval (in Figure 3.1 
.y = a +J ■ is the intersection of the stripe defined by the interval evaluations of 
f(a) and f(b) with the real axis). If the width w(J) of J is more than 4, we set p to 
2p and retry. Otherwise, let £ be the integer closest to mid(7) and set m* :=a + £- 

For m = a+ f(a~yj(b) ~ fl ) as ^ )er ' ore an< i m j := a + J ' >2 W *- rec ' dots) for j = 0, . . . ,N, 
the following Lemma shows that the computed m* = mi indeed approximates m on the 
m ; -grid: 

Lemma 1. Let m be inside the subinterval \mj,trij+\\. Then, m* = mj or m* = 
Moreover, let m! £ {m.j,mj + \\ be the point that is closer to m. If \m — m'\ < then 
m* = m' . 

Proof. Let A := ^ f(a)-f(b) an< ^ ^ t ^ le m t erva l computed by interval arithmetic as above, 
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with width at most \. Since m — f(a) + X^jf- <G [mj,mj + i], it follows that j < X < 



j+l. By construction, X 6 /. Therefore, |A — mid(7)| < g and, thus, it follows that 
mid(7) can only be rounded to /' or j + l. Furthermore, for m' = nij, \m — m'\ < >2 ^- 
implies that |A — j\ < \. It follows that |mid(7) — j'| < I by triangle inequality, so 
mid(y) must be rounded to j. The case m' = nij + \ is analogous. □ 

The second substep to replace in the QIR method is to check for sign changes 
in subintervals in Steps 6|8 As before, we set CO := w(I)/N. Instead of comparing 



the signs at m' and tn! ± CO, we choose the seven subdivision points (red crosses in 



Figure 3.1 



m* — CO,m* — , m* — — , m* , m* + — , m* — , m* + CO. (3.1) 

8 2 2 8 



In case that m* = a or m* = b, we only choose the 4 points of (3.1 1 that lie in /. For a 
working precision p, we evaluate the sign of / at all subdivision points using interval 
arithmetic. If the sign remains unknown for more than one point, we set p to 2p and 
retry. After the sign is determined for all except one of the points, we look for a sign 
change in the sequence. If such a sign change occurs, we set the corresponding interval 
/* as isolating and call the Aqir step successful. Otherwise, we call the step failing and 
keep the old isolating interval. As in the exact case, we square up after a successful 
step, and reduce it to its square root after a failing step. See Algorithm[3]for a complete 
description. 

Note that, in case of a successful step, the new isolating interval I* satisfies g^vv(/) < 
w(I*) < jjw(I). Also, similar to the bisection method, the function can only be zero 
at one of the chosen subdivision points, and the function is guaranteed to be reason- 
ably large for all but one of them, which leads to a bound on the necessary precision 
(Lemma FT}. The reader might wonder why we have chosen a non-equidistant grid 
involving the subdivision points m* ± | CO. The reason is that these additional points 
allow us to give a success guarantee of the method under certain assumptions in the 
following lemma which is the basis to prove quadratic convergence if the interval is 
smaller than a certain threshold (Section[5^2|. 



Lemma 2. Let I — (a,b) be an isolating interval for some root t, of f, s = sign(/(a)) 
and m as before. If \m— % \ < = |?, then AQ\R(f,I,N,s) succeeds. 

Proof. Let m* be the subdivision point selected by the Aqir method. We assume 
that m* {a,b}; otherwise, a similar (simplified) argument applies. By Lemma [T] 
m e [m* — \ co,m* + and, thus, % e {m* — \co,m* + \co). It follows that the leftmost 



two points of (3.1 1 have a different sign than the rightmost two points of (3.1 1. Since 
the sign of / is evaluated for at least one value on each side, the algorithm detects a 
sign change and, thus, succeeds. □ 



4 Analysis of an AQIR step 

The running time of an Aqir step depends on the maximal precision p needed in the 
two while loops (Step|5] Steps [TT|[T5j ) of Algorithm[3] The termination criterion of both 
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loops is controlled by evaluations of the form 58 (J?, p), where £ is some polynomial 
expression and p is the current working precision. 

We specify recursively what we understand by evaluating £ in precision p with 
interval arithmetic. For that, we define down(x,p) for x E R and p e N to be the 
maximal xq < x such that xq — for some integer k. The same way up(.x,p) is the 
minimal xq > x with xq of the same form. We extend this definition to arithmetic 
expressions by the following rules (we leave out p for brevity): 

down(£i +£2) := down(2?i) + down(£2) 
up(Ei+E 2 ) := up(£i)+up(£ 2 ) 
down(£i -£2) := down(min{down(£i)down(£2),up(£i)up(£2), 
up (Ei ) down (£2 ) , down [E\ ) up (£2 ) } ) 
up(£i-£2) := up(max{down(£i)down(£2),down(£i)up(£2), 
up (£1 ) down (£ 2 ) , up (£1 )up (£ 2 ) } ) 
down(l/£i) := down (1 /up (£1)) 
up(l/£j) := up(l/down(£i)) 

Finally, we define the interval 58(£,p) := [down(£,p),up(£,p)]. By definition, the 
exact value of £ is guaranteed to be contained in 58(£,p). We assume that polynomi- 
als / e R[x] are evaluated according to the Horner scheme, and when evaluating f{c) 
with precision p, the above rules apply in each arithmetic step. The next lemma pro- 
vides a worst case bound on the size of the resulting interval 58 (/(c), p) under certain 
conditions. We further remark that, in an actual implementation, 58(£,p) is usually 
much smaller than the worst case bound derived here. Nevertheless, our complexity 
analysis is based on the latter bound. Throughout the following considerations , F G N 
denotes an integer upper bound on the root bound Ff, that is, F>Ff, and, in particular 
log |zj| < r for all roots Zi of /. 



Lemma 3. Let f be a polynomial as in (O), eel with \c\< 2 r+1 , and p e K Then, 



\f{c)-dovm(f(c),p)\<2-P + \d+\) 2 2 x+di - T+2) (4.1) 
|/( C )-up(/(c),p)| <2-"+V+l) 2 2 T+</(r+2) (4.2) 

In particular, 58 (/(c), p) has a width of at most 2~ p+2 (d + i~) 2 2 z+d ( r+2 \ 

Proof. We do induction on d. The statement is clearly true for d = 0. For d > 0, we 
write /(c) = uq + cg(c) with ao£M the constant coefficient of / and g of degree d—1. 
Note that, for any real value x, |down(x,g) — x\ < 2~ p , same for up. Therefore, we can 
bound as follows (again, leaving p out for simplicity): 

|/(c)-down(/(c))| = \a +cg(c)-down(a + cg(c))\ = \a + cg(c) -down(a ) -down(cg(c))| 
< |cg(c)-down(cg(c))|+2~ p 

Note that down(c • g(c)) = down(//i (c) • H2(g(c))) where H1.2 = down or H\ 2 = up. 
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Algorithm 3 Approximate Quadratic interval refinement 



Input: / e M.[x] square-free, / = (a,b) isolating, JV = 2 2 'eN,s = sign {f(a)) 
Output: (J,N r ) with J CI isolating and iV'eN 
l: procedure Aqir(/,7 = (a,b),N) 

2: if N = 2, return (Approximate_Bisection(/,/,.s),4). 
3: CO <- ^ 
4: p 4-2 

5: while 7 4- <B(N f ^"j^ ,p) has width > i, set p 4- 2p 
6: m* 4- a + round(mid(7)) • 0) 

7: if m* = a, s 4- 4,V 4- [m*,m* + j(B,m* + |a>,w* + Co],S 4- [s, 0,0,0] 
8: if m* =b,s(- 4,V 4- [m* - 0),m* - |t»,w* - \(Q,m*],S 4- [0,0,0,— s] 
9: if a < w* < b, s 4— 7,V 4- [m* — CO,m* — |ei),m* — ^CO,m*,m* + \co,m* + 
l(0,m* + co],S^- [0,0,0,0,0,0,0] 

10: p4-2 

11: while S contains more than one zero do 
12: for i=l,. . . ,s do 

13: If S\i] = 0, set S\i] 4- sign <B (/( V [i] ) , p ) 

14: end for 

15: P^-2p 

16: end while 

17: If 3v,w : 5[v] -S[w] = -1 A(v+1 = ivV (v + 2 = wAS[v+l] = 0)) return 

((V[v],V[w]),JV 2 ) 
18: Otherwise, return (/, v^V) 
19: end procedure 



Moreover, we can write H\(c) = c — e with |e| < 2 p . Therefore, we can rearrange 

\cg(c) -down(eg(c))| +2~P < \cg(c) - (c - e) ■ H 2 (g(c))\+2-P +1 

< \cg(c) - c ■ H 2 (g(c))\ + \e\- \H 2 (g(c))\ + 2-P+ 1 

< \c\ ■ \g(c)~H 2 (g(c))\+2-P\H 2 (g(c))\+2-P+ 1 

By a simple inductive proof on the degree, we can show that both |up(g(c))| and 
|down(g(c))| are bounded by d2 r+d ( r+2 \ Using that and the induction hypothesis 
yields 

\c\-\g(c)-h(g(c))\+2-P\H 2 (g(c))\+2-P +1 

< 2 r+2 2~ p+1 t/ 2 2 T+(rf ~ 1)(r+2) + 2~ p d2 T+d( - r+2) +2~ p+1 

< 2-P +l {d 2 + d+\)2 z+d V + V <2-P + \d+\) 2 2 xd 

The bound for |/(c) — up(/(c)) | follows in the same way. □ 

For the sake of simplicity, we decided to assume fixed-point arithmetic, that means, 
p determines the number of bits after the binary point. We refer the interested reader 
to lf]~3l Thm. 12], where a corresponding result for floating-point arithmetic is given. 
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We analyze the required working precision of approximate bisection and of an 
Aqir step next. We exploit that, whenever we evaluate / at t subdivision points, t — 1 
of them have a certain minimal distance to the root in the isolating interval. The follow- 
ing lemma gives a lower bound on |/(*o) I f° r sucn a point xq, given that it is sufficiently 
far away from any other root of /. 



Lemma 4. Let f be as in (1.1 \, 4 = z, a real root of f and xq be a real value with 
distance \xq — Zi\ > -/ to all real roots Zi ^ Zi - Then, 

\f(x )\>\^-x \.2-^^f). 
( recall the notations from Section^for the definitions of Oi and T.f) 

Proof. For each non-real root zi of /, there exists a complex conjugate root z< and, thus, 
we have |jco — Zi\ > Im(z,) > y > ^ for all i = m + 1, . . . , d as well. It follows that 

d 

|/(*o)| = \a d Y\{xQ-Zi)\ = \a d \ ■ |4 -xq\ ■ \\ \x -Zi\ 

i=l i=l,... : d:i^i(j 

>\t-M-^--n^>\?;-M-2- 2d - r -2- z f, 

a 'o i=l ^ 



where the last inequality uses that |z,-| < 2 r and, thus, o"(z,) < 2 r+1 . □ 
We next analyze an approximate bisection step. 



Lemma 5. Let f be a polynomial as in (1.1 1, / = (a,b) C (— 2 r+2 ,2 r+2 ) be an isolating 
interval for a root 4 = z, off and s = sign(/(a) ). Then, Algorithm^applied on (/, /, s) 
requires a maximal precision of 

p :=2\og{b-a)- x + 41og(d + l) + 8d + 10+2(d+l)r+T + 2£ / 
= 0(log(fc - a)' 1 + T + dT + £/), 

anj ;fs fe/f complexity is bounded by 0(d(\og(b — a) -1 + T + rfr + Zy)). 

Proof. Consider the three subdivision points mj :=«+/'• ^p, where 1 < j < 3, and 
an arbitrary real root z ( - ^ ^ of f. Note that |m 7 - — z, | > because the segment from 
m ; - to Zi spans at least over a quarter of Moreover, |4 — m ; | < — a), and so 

3 

O; < 14 -z/l < + < j(b-a) + \mj-Zi\<4\mj-zi\. 

It follows that mj has a distance to zi of at least ^ . Hence, we can apply Lemma |4] 
to each mj, that is, we have |/(m 7 )| > |4 — my| ■ 2~( 2i/+r+J -/). Since the signs of / at 
the endpoints of / are known, it suffices to compute the signs of / at two of the three 
subdivision points. For at least two of these points, the distance of mj to 4 is at least 
^p, thus, we have \f{mj)\ > \b -a\ ■ 2~( 2d+3+r+ ' L f) for at least two points. Then, due 
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to Lemma [3] we can use interval arithmetic with a precision p to compute these signs 
if p satisfies 

2-" +2 ( d /+l) 2 2 T+ ' /(r+2) < (b-a)-2-( 2d+3+r+z f\ 

which is equivalent to p > ^. Since we double the precision in each step, we will 
eventually succeed with a precision smaller than po. The bit complexity for an arith- 
metic operation with fixed precision p is 0(p + dz). Namely, since the absolute value 
of each subdivision point is bounded by 0(t), the results in the intermediate steps 
have magnitude 0(dt) and we consider p bits after the binary point. At each sub- 
division point, we have to perform 0(d) arithmetic operations for the computation 
of f(ntj), thus, the costs for these evaluations are bounded by 0(d(dz + p)) bit op- 
erations. Since we double the precision in each iteration, the total costs are domi- 
nated by the last successful evaluation and, thus, we have to perform O(d(po +dx)) — 
d(d(log(b-a)~ l +dT + T.f)) bit operations. □ 

We proceed with the analysis of an Aqir step. In order to bound the required 
precision, we need additional properties of the isolating interval. 



Definition 6. Let f be as in {1.1 1, / := (a,b) be an isolating interval of a root £ of f. 
We call I normaQ;/ 

. /C(-2 r+2 ,2 r + 2 ), 

• \P ~~ z i I > T f or evel J P ^ I an d Zi 7^ <=> an d 

. min{|/(a)|,|/(/7)|} >2-( 28+2T +™ r+2I /- 51o ^- fl )). 

In simple words, a normal isolating interval has a reasonable distance to any other 
root of /, and the function value at the endpoints is reasonably large. We will later see 
that it is possible to get normal intervals by a sequence of approximate bisection steps. 

Lemma 7. Let f be a polynomial as in \1-1\ , I — (a,b) be a normal isolating interval 
for a root £ = z k off with s = sign(/(a)), and let N < 2 2 ( r + 4 - lo s( i - a )). Then, the 
AQIR step for (f,I,N,s) requires a precision of at most 

Pmax ■= 87c/T+ HdT + AZf- 141og(fc- a) 

and, therefore, its bit complexity is bounded by 

d(d(T + dr + Z f -log(b-a))). 

Moreover, the returned interval is again normal. 

Proof. We have to distinguish two cases. For > 2, we consider the two while-loops 
in AlgorithmjiJ In the first loop (Stepjijl, we evaluate jV f^a)-\(b) v ^ a mterva l arithmetic, 

1 The reader may notice that the definition of "normal" depends on the upper bound T on Vf . Throughout 
our argument, we assume that such an initial T is given. We will finally choose a T which approximates Tf 
up to an (addative) error of O(logrf). 
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doubling the precision p until the width of the resulting interval J is less than or equal 
to 1 /4. The following considerations show that we can achieve this if p fulfills 

2-P+ 2 {d + 1 ) 2 2 T +<'( r + 2 > < min (l/( fl )U/(fr)l) _ (4 3) 

W.l.o.g., we assume f(a) > 0. If p fulfills the above condition, then, due to Lemma[3] 
Q3(iV-/(a),p) is contained within the interval 

[Nf{a) _ \lM iNf(a) + =Nf[a) . (1 _ J_, ! + J_] 

32 32 32A? 32JV J 

and 03 (/(a) —f(b), p) is contained within the interval 

[/(«) -/(*) - ,/(«) - /(*) + - (/(«) -/(*)) ■ [1 - 3^,1 + 3^], 

where the latter result uses the fact that /(a) and f(b) have different signs. It fol- 
lows that Q3(AT /(a ffi (fe) i p) is contained within • [(1 - g^)/(l + g^), (1 + 
32jv)/(l — 52v)]' an d a simple computation shows that N- [(1 — 32v)/(l + 557)1 (1 + 
3SvV — 3Sv)] ^ as width less than 1/4. Hence, since f(a) fib) ^ as absolute value less 
than 1, 03 
writes as 



than 1, ^(^ /(a)-/(A) >P) nas width less than 1 /4 as well. The bound (4.3 1 on p also 



p>l + 2\og(d+\) + x + dT + 2d + \ogN + \ogrmn(\f{a)J{b)\)- 1 

■ M 

M-f{b) 



and since we double p in each iteration, computing N via interval arithmetic 



up to an error of 1/4 demands for a precision 

p < 14 + 41og(t/+l)+2T + 2t/r + 4d/ + 21ogA' + 21ogmin(|/(fl),/(Z7)|)- 1 
< 14 + 2T+10dr + 21ogAr + 21ogmin(|/(a),/(Z?)|)- 1 , 

Since / is normal and because of the posed condition on N, we can bound this by 

p < HdT+4(T+5-log(fe-a))+2(32rfT+2E / -51og(i-fl)) 
< 87dx +4Lf - 141og(fo -a) < p max . 

We turn to the second while loop of Algorithm [3] (Steps [TTfl"5j ) wher e / is evaluated 
at the subdivision points m* — co ,m* — ^ , . . . ,m* + CO as defined in (3.1 1. Since the 
interval is normal, we can apply Lemma [4] to each of the seven subdivision points. 
Furthermore, at least six of these points have distance > to the root B, and, thus, 
for these points, |/| is larger than ■ 2~( M + T+I V). Then, according to Lemma 4.3 it 
suffices to use a precision p that fulfills 

2-" + V+l) 2 2 T+rf(r+2) < ^■^ {74+V+Zf \ or 
P >Pi :=6 + 21og(d+l) + T + c/r + 4o' + r + E / + logjV-log(Z?-a). 
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The same argumentation as above then shows that the point evaluation will be per- 
formed with a maximal precision of less than 

2pi <2(6 + T + 7rfT + r+£/ + logiV-log(£-fl)) 

< 12 + 2T+14ir + 2r + 2£ / + 4(r+4-log(fo-a))-log(&-a) 
<28+2T+17 C /r + 2E / -51og(fe-a) 

which is bounded by p max . Moreover, at the new endpoints a 1 and b ', \ f\ is at least 

2~2pi > 2-(28+2T+17(/r+2I f -51og(6-a)) > 2-(28+2T+17dr+2Z / -51og(6'-a')) 

which proves that /' = (a' ,b') is again normal. 

It remains the case of /V = 2, where a bisection step is performed. It is straight- 
forward to see with Lemma[5]that the required precision is bounded by p max , and in an 
analogue way as for the point evaluations for TV > 2, we can see that the resulting inter- 
val is again normal. By the same argument as in Lemma[5] the overall bit complexity 
of the Aqir step is bounded by 

d(dp max ) = d(d(d'Z+Lf-log(b-a))). □ 



5 Root refinement 

We next analyze the complexity of our original problem: Given a polynomial / as in 



( 1. 1 1 and isolating intervals for all its real roots, refine the intervals to a size of at most 
2~ L . Our refinement method consists of two steps. First, we turn the isolating inter- 
vals into normal intervals by applying bisections repeatedly. Second, we call the AQIR 
method repeatedly on the intervals until each has a width of at most 2~ L . Algorithm[5] 
summarizes our method for root refinement. We remark that depending on the prop- 
erties of the root isolator used to get initial isolating intervals, the normalization can 
be skipped; this is for instance the case when using the isolator from 1171 . We also 
emphasize that the normalization is unnecessary for the correctness of the algorithm; 
its purpose is to prevent the working precision in a single Aqir step of growing too 
high. 



5.1 Normalization 

The normalization (Algorithm |4]i consists of two steps: first, the isolating intervals 
are refined using approximate bisection until the distance between two consecutive 
intervals is at least three times larger than the size of the larger of the two involved 
intervals. This ensures that all points in an isolating interval are reasonably far away 
from any other root of /. In the second step, each interval is enlarged on both sides 
by an interval of at least the same size as itself. This ensures that the endpoints are 
sufficiently far away from any root of / to prove a lower bound of / at the endpoints. 
W.l.o.g., we assume that the input intervals are contained in (— 2 r+1 ,2 r+1 ) because all 
roots are contained in that interval, so the leftmost and rightmost intervals can just be 
cut if necessary. Obviously, the resulting intervals are still isolating and disjoint from 
each other. Moreover, they do not become too small during the bisection process: 
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Algorithm 4 Normalization 

Input: / g R[f] a polynomial as in (JTTTJ, 4 = (ai,Z>i), . . . ,/„, = (a m ,b m ) disjoint iso- 
lating intervals in ascending order, Si,... ,s m with s k = sign(/(min4)) 
Output: normal isolating intervals Ji,...,J m with 6 7^ n 
i: procedure normalize^,/! , . . .,I m ) 
for k=l,. . . ,m-l do 

while min4+i — max4 < 3max{w(4), w(4+i)} do 
if w(4) > w(4 +1 ) 

then Approximate jisection(/, 4,**;) 
else Approximate_bisection(/, 4 +1 ,Sk+i_) 
end while 
end for 

for k=l,. . . ,m-l do 

dk min4 + i — max4 

J k <— [a k — d k _\/3,b k + d k /3] >enlarge4by more than w (4) at both sides 
end for 

return J\,...,J m 
end procedure 



Lemma 8. For J\ , . . . , J m as returned by Alg. Q w(4) > 5 Oft- 

Proof. After the first for-loop, the distance <4 between any two consecutive intervals 
4 and 4+1 fulfills dk > 3max{w(4), w(4+i)}, thus o k < w{I k ) + w(4+i) + d k < 2d k . 
Hence, in the last step, each I k is enlarged by at least (J k /6 on each side. This proves 
that the corresponding enlarged intervals J k have size <J k /3 or more. □ 

Lemma 9. Algorithm^is correct, i.e., returns normal intervals. 

Proof. Let J\,...,J m denote the returned intervals, and fix some interval J k containing 
the root z k of /. We have to prove the three properties of Definition[6] The first property 
is clear because the initial interval are assumed to lie in (— 2 r+1 ,2^ +1 ). In the proof 
of Lemma [I] we have already shown that 4 is eventually enlarged by at least <J k /6 on 
each side. More precisely, the right endpoint of J k has distance at least d k /3 > G k+ i/6 
to J k +\, and the left endpoint of J k has distance at least d k -\/3 > O k -\/6 to J k -\. It 
follows that, for each xq £ J k , we have |xo — z k ±i \ < <J k ±i/3, respectively. Hence, the 
second property in Definition|6]is fulfilled. 

For the third property of Definition [6] let e be one of the endpoints of J k . We have 
just proved that the distance to every root Zi except z k is at least ? and |e — I > Oit/6. 
With an estimation similar as in the proof of Lemma|4] we obtain: 

i/(*)i>fn? >^- 4 i T 2-/=2-(^), 

i^k 

and 2-( M + E / +1 ) > 2-( 28 +2T+™r+2i / -5io g (fo-«)) beca use \og{b-a) < T+2 and -£/ < 
d{V+\)<2dY. □ 



14 



Algorithm 5 Root Refinement 



Input: / = Y, a i x ' G H^[ f ] a polynomial as in (JTTTj), isolating intervals I\, . . . ,I m for the 
real roots of / in ascending order, LeZ 
Output: isolating intervals J\,...,J m with w(J k ) < 2~ L 
i: procedure root_refinement(/,L,4, . . . ,/„,) 

2: s k := sign(a d ) ■ (-l) m - k+1 > s k = sign(/(min/ fc )) 

3: J 1 ,...,J m ^r- NORMALIZE ,...,/„,) 

4: for k=l,. . . ,m do 

5: N ^4 

6: while w(J k ) > 2~ L do (Jk,N) <-AQm(f,J k ,N,Sk) 

i: end for 

8: return J\,...,J m 

9: end procedure 



Lemma 10. Algorithm^has a complexity of 

d(d(dF+Lf)(T;+dF+Lf)) 

Proof As a direct consequence of Lemma [8] each interval 4 is only bisected 0(T + 
log(ajt) -1 ) many times because each starting interval is assumed to be contained in 
(— 2 r+1 ,2 r+1 ). So the total number of bisections adds up to 0(dT + £/•) considering 
all roots of /. Also, the size of the isolating interval 4 is lower bounded by • CTj. = 
2-0(i/+rfr) > so that one approximate bisection step has a complexity of 0(d(x + dT + 
due to Lemma[5] □ 



5.2 The AQIR sequence 

It remains to bound the cost of the calls of Aqir. We mostly follow the argumenta- 
tion from ifTTI . mostly referring to that article for technical proofs. We introduce the 
following convenient notation: 

Definition 11. Let Iq : — I be a normal isolating interval for some real root £ of f, 
No '■= 4 and s := sign(minTo)- The AQIR sequence (Sq,S\ ,...,S v .) is defined by 

So := (7 ,iVo) = (/,4) S, - (7,-,M) := AQi R (/,4_i,iV/_i,i)/or i > 1, 

where is the first index such that the interval I v . has width at most 2~ L . We say that 

Aqir Aqir 
Si — > Si + \ succeeds if AQlR(f,Ii, Ni,s) succeeds, and that Si — > S/+i fails otherwise. 

As in 1111 . we divide the QIR sequence into two parts according to the following 
definition. 

Definition 12. For t, a root of f, we define 



c, 



^l/'^)l+It 2 (^)" 2 l/«(^)l 
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For (So, ■ ■ ■ ,S V . ) the QIR sequence of t,, define k as the minimal index such that S^ = 

(Ik,Nk) A ^> R Sk+i succeeds and w(Ik) < Qj. We call (So, ■ ■ ■ ,Sk) linear sequence and 
(St, . . . , S v . ) quadratic sequence ofE, 

Note that ifTTI defined a different threshold for splitting the qir sequence, and the 
linear sequence was called initial sequence therein. We renamed it to avoid confusion 
with the initial normalization phase in our variant. 



Quadratic convergence. We start by justifying the name "quadratic sequence". 
Indeed, it turns out that all but one Aqir step in the quadratic sequence are successful, 
hence, N is squared in (almost) every step and therefore, the refinement factor of the 
interval is doubled in (almost) every step. We first prove two important properties of 
Qj as defined in Defition 
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Lemma 13. Let B, e C be a root of f. 

1. 0<C^ < -gp- 

2. Let /J. S C be such that |£ — /i| < Q=. Then 

c c 1/(4)1 

1 8|/"(M)I' 



Proof. Note that all summands in the denominator of Q- are non-negative. Therefore, 
the first property follows immediately by removing all but the first summand in the 
denominator. 

For the second property, we consider the Taylor expansion of /"(/i) in % : 
/''( M ) = £( M _|)«^|I. 

Because |ju — £ 1 < Q= < by the first property, we can bound 

\f"w\<i(^y 2 \f {,) m 



It follows that 



v«i 8 ( Et2 (^)'-V>«)i ' 



□ 



The following bound follows from considering the Taylor expansion of / at £ in 
the expression for m: 
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Lemma 14. Bill Thm. 4.8] Let (a,b) be isolating for % with width 8 < Q= and m as 
in Lemma^^f i.e., m = a+ j{Jy~f(g) (b ~ a ) )■ Then, \m — ^ | < 

Proof. We consider the Taylor expansion of / at ^ . For a given x £ (a,b),we have 

f( x )=f'(t;)( x -t;)+ l 2 f"(l)( x -t;) 2 
with some B, £ [x,^] or ^ £ Thus, we can simplify 



f(b)(a-S)-f(a){b-S) 



fQ>)-f(a) 



\f{^){b-^-{a-^)-f"{^){a-^{b-^)) 



f{b)-f{c 



< ^1*— €1 m ^ Ml < WM{ 

for some v £ (a,b). The Taylor expansion of / yields f(v) = f'(^) +/"(v)(v - 1) 
with v G (a,/?). Since 8 < it follows with Lemma 
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|/"(v)(v-§)|<|/"(v)|Q<i|/'(§)|. 



Therefore |/'(v)| > g|/'(§)| > ||/(|)|, and it follows again with Lemma 



13 that 



1 51 - 1/^)1 



g 2 



8 2 

<8c; D 



8max{|/"(^)M/"(fe)|} 



Corollary 15. Let Ij be an isolating interval for c, of width 8j < Then, each call 
of the AQIR sequence 



(Ij,Nj) A T(I j+1 ,Nj +1 )^ 



succeeds. 



Proof. We use induction on i. Assume that the first i Aqir calls succeed. Then, another 
simple induction shows that 5 ;+I - := w(/ ; +,) < ^7 < j^-, where we use that Nj+i = 



Nj +i _ v Then, according to Lemma 
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Nj+i " N j+i > 
we have that 



C{ 1 _ 1 S j+1 



'8Q "= J+r N J+i 8C^ SNj+i 1 
with m as above. By Lemma|2] the Aqir call succeeds. 

Corollary 16. iffi! Cor. 4.10] In the quadratic sequence, there is at most one failing 
Aqir call. 



□ 
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Proof. Let (7;,A/;) A ^> R (/;+i,A/; + i) be the first failing Aqir call in the quadratic se- 
quence. Since the quadratic sequence starts with a successful Aqir call, the predeces- 

It ., \ AQI. 

sor (Ii-i,Ni-x) -> 
have the sequence 



sor (/j-i, iV/-l) A ^> R (Ii>Ni) is also part of quadratic sequence, and succeeds. Thus we 



Sucess Fail 

(I.-uN.-i) A -T (I„N,) A ^ R (I i+ uN i+l ) 



One observes easily that w(I i+l ) = w(/,) = ^ < andA/, + i = ^/N) = yNf_ l = 
Ni-i. By Corollary [T3] all further Aqir calls succeed. □ 

Cost of the linear sequence. We bound the costs of refining the isolating interval 
of % to size Q- with Aqir. We first show that, on average, the Aqir sequence refines 
by a factor two in every second step. This shows in particular that refining using Aqir 
is at most a factor of two worse than refining using approximate bisection. 

Lemma 17. Let (Sq, . . . ,Sj>) denote an arbitrary prefix of the Aqir sequence for 
starting with the isolating interval Iq of width 8. Then, the width ofl? is not larger than 

Proof. Consider a subsequence (5/, . . . ,5,-+ ; -) of (So, ■ ■ • ,Si) such that Si A -^> R S/+i is 
successful, but any other step in the subsequence fails. Because there are j steps in 
total, and thus j —I consecutive failing steps, the successful step must have used a N 
with N > 2 V ~ X . Because V~ l > 5, it holds that 

w(I 1+j ) < ^ < w(I l+J )2- 2J - 1 < W {I i+] )2-Jl\ 

Repeating the argument for maximal subsequences of this form, we get that either 
w{l£) < w{Iq)2~ l I 2 if the sequence starts with a successful step, or w(Ii) < w(Iq)2~^~ 1 ^ 2 
otherwise, because the second step must be successful in this case. □ 

We want to apply Lemma[7]to bound the bit complexity of a single Aqir step. The 
following lemma shows that the condition on N from Lemma [7] is always met in the 
Aqir sequence. 

Lemma 18. Let (Ij,Nj) A ^> R (Ij+i,Nj+i) be a call in an Aqir sequence and Ij :— 
(a,b). Then, Nj < 2 2 ( r + 4 - lo s( fc -")). 

Proof. We do induction on j. Notethat/oC (— 2 r+2 ,2 r+2 ) by normality, hence b — a < 
2 r+3 . It follows that 2 2 ( r+4 - lo §( /: '- a )) > 4 = Nq. Assume that the statement is true for 
j — 1. If the previous step (7,_i,A/,_i) -^> R {Ij,Nj) is failing, then Nj = yjNj- \ and 
the isolating interval remains unchanged, so the statement is trivially correct. If the 
step is successful, then it holds that (b — a) < ^7=- By rearranging terms, we get that 

jy. < 22(r+3-log(/j-a)) q 
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It follows inductively that the conditions of Lemma [7] are met for each call in the 
Aqir sequence because Iq is normal by construction. Therefore, the linear sequence 
for a root % of / is computed with a bit complexity of 



oar+iog^rvaog^+T+dr+E/)) 



(5.1) 



because (9(r + log(C^ 1 )) steps are necessary to refine the interval to a size smaller than 
Qj by Lemma|l7 and the bit complexity is bounded by (9(c/(log(C^ 1 ) + X + dT+'Lf)) 
with Lemma l7llt remains to bound log(C,c) _1 ; we do so by bounding the sum of all 
log(C,c) _1 with the following lemma. 

Lemma 19. ZT= l lo g(Q,) _1 = 0(d{F + logd) +£/)) 
Proof. We note that 











\°e f 






)) 



We focus on the quotient 



Let Zj, . . . ,z' d _ l denote the (not necessarily distinct) 



roots of /'. Note that for ieC and any i > 1, 

f { Hx)=a d £ J! 



xc{l n-1] je{l d-1) 

\X\=i-l j^X 



Therefore, the quotient writes as 

f {,) (zi) 



f'(ze) 



e n 



xc{i,...,tf-i} jex z t z j 
pt|=i-l 



xc{i,....rf-i} yeX I** 
1X1=1-1 



Since — > ^ Thm.8], we can further bound this to 



e n^< e n^< e 

xc{i,...,i-i} ;eX I** xc{i i-i> jex xc{i,...,tf-i} 

ixi=i-i ix"i=f— i |xi=;-i 



i-l 



and, therefore, 



d 2 J 



< 



d 'a e y-2d 2i - 2 

r i-l 



d 



d J 2 



i-l 



i-l 



= E^ = (^-D 



Plugging in into the overall sum yields 

m 



" / /rf 2 rf 2 



m Si 

= 3</+£log- 



= 3d + 3mlogd + Zf+ E lo g°f < 3d + 3dlogd + T. f + d(T + 1) = 0(<i(r+logd) +£/)□ 

£=m+l 
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Lemma 20. The linear sequences for all real roots are computed within a total bit 
complexity of 

6{d{dY+Z f ){x + dY+Z f ) 
Proof. The total cost of all linear sequences is bounded by 

m 

6(£ (r + io g (c- 1 ))d (log^ 1 ) + t + dY + z f )) . 

By rearranging terms, we obtain 

= d{d 2 r(T+dr+i: f )+d(T+dr+z f )^iog(c- 1 )+d{Y^iog(c- 1 ))^ 

which equals O (d (dY + Zf) (t + dY + Zf) ) with Lemma 
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□ 



Cost of the quadratic sequence. Let us fix some root | of /. Its quadratic se- 
quence consists of at most 1 + logL steps, because N is squared in every step (except 
for at most one failing step) and the sequence stops as soon as the interval is smaller 
than 2~ L . Since we ignore logarithmic factors, it is enough to bound the costs of one 
QIR step in the sequence. Clearly, since the interval is not smaller than 2~ L in such a 
step, we have that \og(b — a) < L. Therefore, the required precision is bounded by 
0(L + x + dT + Zf) . It follows that an AQIR step performs up to 0(d(L + T + dY + Zf)) 
bit operations. 

Lemma 21. The quadratic sequences for one real root is computed within a bit com- 
plexity of 

6(d(L + T + dY+Z f )). 
Total cost. We have everything together to prove the main result 
Theorem 22. Algorithm^performs root refinement within 

d(d(dr f + Z f ) 2 +dL) 
bit operations for a single real root^of f, and within 

d(d(dT f + Z f ) 2 +d 2 L) 

for all real roots. The coefficients of f need to be approximated to 0(L + dTf + Zf) 
bits after the binary point. 

Proof. We first restrict to the case where 1 < < 2. The so far achieved complex- 
ity bounds are formulated in terms of an arbitrary (but given) upper bound Y € N on 
ry. In (T7, Section 6.1], it is shown how to compute a Y with Yf < Y < Yf + 41ogaf 
using 0((dYf) 2 ) bit operations and approximations of / to O(dYf) bits after the bi- 
nary point. Furthermore, the latter construction also shows that T = [log(max,- |a;Q] = 

2 In its initial formulation, Algorithm[5]assumes that isolating intervals for all real roots are given. If only 
one isolating interval 4 for a root Zk is given, we have to normalize 4 first and, then, compute the signs of / 
at the endpoints of /. 
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0(dT) if 1 < \a c [\ < 2. By Lemma [Toj the normalization for all isolating intervals 
requires 0(d(dT + Hf)(t + dT + Hf)) bit operations. The linear subsequences of the 
Aqir sequence are computed in the same time by Lemma |20| The quadratic subse 
quences are computed with 0(d 2 L + d 2 x + d 3 Y + d 2 T.f) bit operations by Lemma 21 



the latter three terms are all dominated by 0(d(dY + "LA{x + dT +£/)). Hence, with 
r = 0{Tf + \ogd) as above and T = 0(dYf), the claimed bound on the bit complex- 
ity to refine all roots follows. The maximal number of required bits follows from 
Lemma [7] because the maximal required precision in any AQIR step is bounded by 
0(L + 1 + dF + £/) = 0{L + dTf +£/). The bound on refining a single root follows 
easily when considering the cost of the quadratic sequence for this root only. 

For the more general case, where 1 < \a c [\ < 2 is not necessarily given, we first 
shift the coefficients by s = [log \a^\\ bits such that we can apply the above result to 
the shifted polynomial. Since this coefficient shift does not change the roots, our bit 
complexity bound follows immediately. For the required precision, we need 0(L + 
dFf + Lf) — s since we need an approximation of the shifted polynomial to 0(L + 
dTf+Y,f) bits after the binary point. □ 

For integer polynomials, we have Ff = 0(t) and = 0(dz) [17, §7.2]. Thus, it 
follows 

Corollary 23. If f is a polynomial with integer coefficients of maximal bitsize X, the 
bit complexity of Algorithm^is bounded by 

0(dh 2 +d 2 L). 

This improves the bound from [11] by a factor of d. 



6 Concluding Remarks 

We have presented a complete solution to the root refinement problem using validated 
numerical methods in this paper. Despite the relative simplicity of the approach, we 
obtain a bit complexity which is essentially competitive to best known bounds which 
have been achieved by much more sophisticated algorithms. 

We have shown that the complexity of approximating roots of a real polynomial 
only depends on the geometry of the roots and not on the complexity or the type of the 
coefficients. By means of the following example, we demonstrate how to benefit from 
such an approach. Assume that the root refinement problem is applied to a non-square- 
free integer polynomial F of degree d and coefficient size T. It is known that its square- 
free part / can be computed in 0{d i x) operations using the Euclidean algorithm and / 
is of degree at most d and has coefficient size at most T + d. However, the roots of F 
and / coincide, so it is possible to apply Cauchy's bound on F (instead of on f) which 
yields a root bound Y € 0{%) (in comparison to Y E 0(x + d) when Cauchy's Bound 
is applied to the square-free part). This finally leads to a complexity of d{d i T 2 +d 2 L) 
for non-square-free polynomials as well. 

Although the focus of this work was the asymptotic complexity, the presented algo- 
rithm also aims for a practically efficient solution of the root approximation problem. 
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Indeed, a simplified version of our approach (for integer coefficients) is included in 
the recently introduced CGAlJ^J-package on algebraic computations 0. Experimental 
comparisons in the context of |2| have shown that the approximate version of QIR 
gives significantly better running times than its exact counterpart. These observations 
underline the practical relevance of our approximate version and suggest a practical 
comparison with state-of-the-art solvers as further work. 
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